Lifetime enhanced transport in silicon due to spin and valley blockade 
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We report the observation of Lifetime Enhanced Transport (LET) based on perpendicular valleys 
in silicon by transport spectroscopy measurements of a two-electron system in a silicon transistor. 
The LET is manifested as a peculiar current step in the stability diagram due to a forbidden 
transition between an excited state and any of the lower energy states due perpendicular valley 
(and spin) configurations, offering an additional current path. By employing a detailed temperature 
dependence study in combination with a rate equation model, we estimate the lifetime of this 
particular state to exceed 48 ns. The two-electron spin- valley configurations of all relevant confined 
quantum states in our device were obtained by a large-scale atomistic tight-binding simulation. 
The LET acts as a signature of the complicated valley physics in silicon; a feature that becomes 
increasingly important in silicon quantum devices. 

PACS numbers: 71.55.Cn, 03.67.Lx, 85.35.Gv, 71.70.Ej 



Silicon, the most important material in the electron- 
ics industry till date, is starting to make its mark in 
nanoscale quantum electronics. While silicon nanowires 
are showing promise in thermoelectric applications [1], 
silicon based qubits [2H6] have been widely investigated 
worldwide due to their long spin coherence times [7] and 
compatibility with the present-day infrastructure of the 
semiconductor industry. In recent times, both top-down 
and bottom-up approaches have been successful in build- 
ing single to few electron silicon devices [8UTU]. The 
demonstration of a single-shot readout of an electron 
spin [11] marks a crucial step towards building a sili- 
con qubit. Unlike most materials used for qubits, silicon 
posses a complex momentum space structure, primarily 
arising due to the six degenerate conduction band (CB) 
valleys away from the Brillouin Zone center. The CB 
valley degeneracy provides an extra degree of freedom to 
the quantum confined orbital states in silicon, and gives 
rise to a variety of novel phenomena that strongly affect 
the electronic structure, transport, and relaxation mech- 
anisms [T2UT6] . It is therefore important to understand 
and quantify valley physics and its effect on electronic 
states. 

In this letter, we report observations of a novel val- 
ley induced phenomena in electronic transport in silicon. 
Through transport spectroscopy measurements of donor 
states in FinFETs, we show that under certain condi- 
tions relaxation of excited states into lower manifolds is 
suppressed due to a combination of both spin and valley 
blockade. This enhanced lifetime results in an additional 
transport path through the excited state, and appears 



as a current step in the stability diagram. The phenom- 
ena dubbed as 'lifetime enhanced transport' (LET) was 
first observed in a silicon double quantum dot [10] due 
to a blocked relaxation of a spin triplet into a ground 
state spin singlet, arising from the long spin relaxation 
times in silicon. In this experiment, LET enables us to 
identify a blocked transition between states that have 
different valley symmetries. We confirm this observa- 
tion (i) by extracting the tunnel rates in and out of the 
donor states through a temperature dependent measure- 
ment and analysis, and (ii) by computing the low-energy 
two-electron spectrum of the system from a multimillion 
atom tight-binding method to compare and identify the 
measured excited manifolds. 




FIG. 1: a) Scanning electron micrograph (SEM) image of the 
FinFET device used in the experiment, b) A ID Schematic 
of the donor potential under an applied electric field. 

The devices used in this study are Si FinFET tran- 
sistors, consisting of a Si nanowire with a gate covering 
three faces of the body (Fig. la). We have measured the 
low-temperature electronic transport of devices where 
the sub-threshold characteristics are dominated by a sin- 
gle As donor atom in the channel. The characteristic fin- 



2 



gerprint of a single dopant consist of a pair of resonances 
associated with the one electron neutral donor (D°) and 
two-electron negatively charged donor (D~) states, a D° 
binding energy of ~ 50 meV, a D" charging energy be- 
tween 30-35 meV, and a spin odd-even effect [19 j. The 
high electric field transforms the confining potential to 
a (hybridized) mix between the donor's Coulomb poten- 
tial and a triangular well at the interface, as depicted 
in Fig. lb) [T7J [18]. This system is essentially a gated 
donor where the donor-bound electrons are partly pulled 
toward the Si/Si02 interface. LET is observed in the 
transport through the two-electron D~ spectrum of the 
donor, as we will show below. 

In a previous work [20 , we measured the Stark shifted 
D° excited state spectrum of an As donor in a number 
of devices, and through atomistic device modeling we in- 
ferred the precise locations of the donors and the fields 
they experienced. Although the charging energy (CE) of 
a bulk As donor is about 52 meV (corresponding to a 
binding energy of 2 mev for D~), the CEs in these de- 
vices are reduced to 30-35 meV due to applied fields and 
capacitive coupling to the gates [19]. As a result, excited 
D~ manifolds form below the conduction band, and are 
visible in the stability diagram. 

Fig. [2^i and(2}} show the stability diagram around the 
D~ and D° charge state respectively. We observe several 
steps in the conducting region which can be associated 
with the excited state spectrum [20j [23], indicated by 
the green dashed lines. These steps are clearly present 
for Vsd < but are absent for Vsd > (except for one 
current step indicated by the red dashed line in Fig. [2^i, 
an indication of strong asymmetric coupling with respect 
to the source and drain contacts . The peculiar current 
step in Fig. [2^i marked by the red dashed line for Vsd > 
is the feature we will focus on. The only way it can 
occur is if the excited state has an unusually long life- 
time so that it does not relax to the ground state during 
transport (Fig. [2]i) [23]. 

The position of the steps in the conducting regions of 
the Coulomb diamond reflect the energy of the excited 
states. The Vsd < trace of Fig. ^ shows that there 
is an excited state between 7.5 meV and 15 meV and an- 
other one at around 20 meV. However, such steps do not 
necessarily have to originate from a single state, a set 
of quantum states can be too close in energy space to be 
distinguished separately. Actually, the second step in the 
conducting region is very broad, starting at 7.5 meV and 
reaching its maximum height at around 15 meV (Fig. 2c), 
which points to an ensemble of states entering the bias 
window. This is confirmed by the trace at Vsd > where 
the single step is much narrower than at the other po- 
larity and must belong to one of the states of this larger 
ensemble. As we will show by an analysis of the tem- 
perature dependence, the first step associated with the 
ground state also originates from multiple states; there is 
at least one first excited state at ~1 meV. These experi- 




FIG. 2: Stability diagrams of the a) D~ and b) D° charge 
state resonance. The green dashed lines indicate current steps 
associated with an excited state entering the bias window. 
The highly asymmetric coupling to the leads and the sequen- 
tial nature of the transport quench the current steps in the 
Vsd > conducting region. The current step in the Vsd > 
indicated by the red dashed line in a) must have a life-time 
(t) > 48 ns to show up significantly in the DC transport, c) 
Source/drain current through the D~ versus energy E (with 
respect to the ground state) for Vg — -30rheV and Vg — 
+30meV, showing the additional current step around 11 meV 
due to LET. d) Schematic of the transport processes through 
the donor atom. Sequential processes are characterized by re- 
laxation to the ground state before tunneling out of the atom, 
indicated by the transport channel r^ 1 . For LET processes, 
such as the one marked by the red dashed in line in a), the 
additional transport path T^ 3 comes into play. 

mental values correspond quite well with the energy spec- 
trum of the two-electron donor calculated from a tight- 
binding based Hartree Fock approach, as shown later. 

A temperature dependent analysis of the ground state 
(GS) conductance peak (Gmax) helps us to verify that 
the transport is in the sequential regime, and enables us 
to extract the tunneling rates Tl and Tr between donor 
and the left and right leads (Fig. |2]i) [21]. Fig. [3}i shows 
the GS resonance at various temperatures. The reso- 
nance broadens with increasing temperature due to the 
temperature-dependent Fermi-Dirac distribution in the 
leads. The resonances were fitted with Gaussian func- 
tions to extract the the peak height (Gmax) and the Full- 
Width-Half-Maximum (FWHM), as shown in Fig.|3]D and 
Fig. [3]3 respectively. Both figures show a transition in the 
transport regime roughly around 5K. Below this transi- 
tion temperature Tc, Gmax increases linearly with 1/T, 
and the FWHM increases linearly with T. This specific 
temperature dependence fully characterizes the transport 
as being sequential and dominated by a single level (i.e. 
the GS). The slope of the FWHM as a function of T is 
equal to 3.2±0.2/c£, very comparable to the theoretical 
value of 3.5 ks [2]. The small discrepancy is attributed 



3 



0. 
0. 

CO 
O 

0. 
0. 
0, 





1.39 K 




1.41 K 




1.5 K 




1.6 K 




1.8 K 


k 


— 2 K 




2.5 K 


o\ 


3 K 




3.5 K 


oi 


4 K 




5 K 




6 K 




7 K 

8 K 




9 K 

10 K 




12.5 K 




15 K 




20 K 




10 15 
T(K) 



FIG. 3: Temperature dependence of the ground state res- 
onance, a) Peak shape at various temperatures, b) Peak 
height (Gmax) and c) FWHM of the resonance as obtained 
by a Gaussian fit. The change in behavior of b) and c) at the 
transition temperature Tc = 5 K shows the crossover from 
single to multi- level transport regime |21j . 

to a finite lifetime broadening in our system. Above the 
transition temperature, the slope of the FWHM as a 
function of T is equal to 4.16 ± 0.02 fc^, comparable to 
the theoretical value of 4.35 ks [H] for transport dom- 
inated by multiple levels. The small discrepancy is at- 
tributed to the fact that the theoretical value holds for 
equidistant level spacing, which is not the case in our 
system. The transition temperature corresponds to the 
first excited state starting to contribute to the transport, 
which yields a splitting between the GS and first excited 
state of ksTc ~ 1 meV. 

We can also obtain the coupling of the GS to the 
source and drain leads from the temperature dependence 
of (Gmax below Tc, which is given by [22| [23] ) 
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By fitting Gmax below Tc to eq. I] we find 
4.2 x 10 13 . Next, we have estimated T(= T L + 
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transport measurements at 300 mK, where T > /c#T and 
the transport is thus in the coherent regime. Here the 
FWHM of the resonance is mainly lifetime-broadened, 
i.e. FWHM = T [25 , and we found T = 173 fieV. Com- 
bining the latter two pieces of information yields Yl = 
173 fieV and Tr = 86 neV. This shows that the donor 
state in this device is coupled asymmetrically with the 
source and drain, with a ratio Tl/Tr w 2000. 

A rate equation model now verifies the observed volt- 
age asymmetry and identifies the peculiar current step of 



Fig. [2^i as arising from LET. In the sequential transport 
regime, the current through the donor atom is given by 



I = e 
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where the subscript denotes the tunneling rate in or out 
of the donor, the superscript indicates the energy level (1 
represents the GS) and n indicates the number of states 
in the bias window defined by the source-drain bias. Thus 
the electron can enter the donor through any number of 
states in the bias window, but leave only through the GS, 
as it eventually relaxes to the GS before exiting. Since 
the donor in our sample is very asymmetrically coupled 
to the leads, r in >> r out for Vsd > while r in << r out 
for Vsd < 0. In this asymmetric limit Eq.|2] reduces to 
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(3) 
(4) 



This shows that if Vsd < 0, the current steps increase as 
more states enter the bias window. On the other hand, if 
^sd > 0, the current becomes saturated and no current 
steps are expected. Actually, with our degree of asymme- 
try the maximum magnitude of any current step (n > 1) 
in Vsd > equals 0.5 % of the step at n — 1, regardless of 
T n . The latter follows directly from Eq.[2]in combination 
with the asymmetry Tl/Tr & 2000 we found from the 
temperature dependent data (Fig. [3|. Thus the current 
step due to LET for Vsd > 0, shown in Fig. [2^i as the 
red line marks a deviation from the sequential transport 
regime, as an extra current path opens up through an 
excited state. By the magnitude of the current steps at 
^sd < we can (conservatively) estimate T^ ~ T^, thus 
the life-time r > 48 ns. 
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FIG. 4: Calculated spin-valley configuration of the D~ level 
spectrum, a) Result of the Hartree-Fock tight-binding cal- 
culations and comparison with the experimental data. The 
energies are labeled relative to the 2e ground state energy, 
b) The spin( up/down) and valley (odd/even) configuration 
of le states showing the generated 2e configurations. From 
D° transport spectroscopy of the system, the donor depth 
and field have been determined to be 4.7 nm and 23 MV/m 
respectively |20| . 
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To identify the spin and valley configurations of the D~ 
excited states, we performed a two-electron (2e) Hartree 
Fock calculation using good one-electron (le) wavefunc- 
tions obtained from atomistic tight-binding simulations 
[26] of the device. In the bulk, a donor ground state has 
equal contributions from all six CB valleys, however, the 
donor in this sample is under a strong z-directed field, 
which makes the k z valley contributions dominant in the 
lowest few le wavefunctions, as we also verified by com- 
puting the Fourier Transform of the wavefunctions. Since 
there is strong valley-orbit interaction arising from the 
sharp interfacial potential and the central-cell donor po- 
tential, the two k z valley states are in symmetric (even) 
and anti-symmetric (odd) superposition [18] . The lowest 
three orbital states shown in Fig. 4a as |0), and 
| A) correspond to even, odd, and even valley symmetries 
respectively. 15 different 2e configurations (i.e. Slater 
Determinants) can be formed from these lowest three or- 
bital (six spin) states, their spin symmetries being either 
singlet or triplet. Some of the non-degenerate 2e config- 
urations are shown in Fig. 4a, marking the occupation 
of the electrons in the spin-valley configurations. 

In Fig. 4b, we compare the 2e spectrum from the- 
ory and measurement. Consistent with the measured 
data of Fig. 2b and 2c, we obtain manifolds of states 
around 1 meV, 8 meV, and 11 meV, each manifold con- 
taining both singlet and triplet states. The lower man- 
ifold of states (< 2meV) does not contain a state with 
a triplet spin- and same (even-even or odd-odd) valley 
-combination, as a triplet spin is not possible within 
the same valley. The first triplet spin- same valley - 
configuration state in the spectrum is theoretically at 
11.2 meV (|A)p) + |D)|A)|T) in Fig. 4b). This particular 
state can not relax back to any of the lower states in the 
spectrum, before either flipping spin (triplet to singlet) 
or flipping valley configuration (even-even to odd-even). 
The | A)|D) + |n)| A)|T) states are very close in energy to 
the experimental lifetime enhanced state at 10 meV. Thus 
the observed LET is due to a state with triplet spin- and 
same valley -configuration, present in the higher mani- 
fold of states. The state is thus life-time enhanced by 
both spin- and valley -blockade. 

In the bulk, the spin relaxation time of donor-bound 
electrons if of the order of milli- seconds, however, the 
vicinity of the Si/Si02 interface is expected to reduce 
that time-scale. We can put a (conservative) lower bound 
on the singlet-triplet relaxation in our D~ system by 
comparing it to the relaxation time Ti of 400 nm quan- 
tum dots defined in a two-dimensional Si/Si02 electron 
system measured by spin echo [27] , which was about 
460 ns. In a silicon 2D confined geometry, valley-lifetimes 
in the order of /is are to be expected [13]. This is all 
consistent with the lower bound of 48 ns valley- lifetime 
estimated here. 

In this letter, we reported the observation of LET in 
silicon arising both due to spin and valley selection rules. 



We showed that in the sequential transport regime, LET 
can be manifested as peculiar current steps in the stabil- 
ity diagram as excited states fail to relax to lower states of 
different valley symmetries, resulting in additional trans- 
port paths. The unusually long-lived states of different 
valley symmetries allude to the fact that the valley index 
can form a two-level system and thus can act as a pseudo- 
spin quantum number. It is interesting to note that this 
pseudo-spin has recently been proposed as the basis in 
a valley-based quantum computation scheme [16] . The 
LET thus serves as a signature of of the complex mo- 
mentum space of silicon, and highlights the importance 
of understanding valley physics for silicon quantum elec- 
tronics. 
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